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Abstract 

Isolated photon production in deep inelastic ep scattering has been measured 
with the ZEUS detector at HERA using an integrated luminosity of 320pb _1 . 
Measurements were made in the isolated-photon transverse-energy and pseudo- 
rapidity ranges 4 < Kf < 15 GeV and —0.7 < rf 1 < 0.9 for exchanged photon 
virtualities, Q 2 , in the range 10 < Q 2 < 350 GeV 2 and for invariant masses of 
the hadronic system Wx > 5 GeV. Differential cross sections are presented for 
inclusive isolated photon production as functions of Q 2 , x, Ej, and rp . Leading- 
logarithm parton-shower Monte Carlo simulations and perturbative QCD pre- 
dictions give a reasonable description of the data over most of the kinematic 
range. 
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1 Introduction 



In the study of high-energy collisions involving hadrons, events in which an isolated high- 
energy photon is observed provide a direct probe of the underlying partonic process, since 
the emission of these photons is unaffected by parton hadronisation. 

Isolated high-energy photon production has been studied in a number of fixed-target and 
hadron-collider experiments [TJ. Previous ZEUS and HI publications have also reported 
the production of isolated photons in photoproduction [2H6] , in which the exchanged pho- 
ton is quasi- real (Q 2 ~ 0), and deep inelastic scattering (DIS) [7J[B], in which Q 2 ps GeV 2 . 

Isolated photons are produced in DIS at lowest order in QCD as shown in Fig. [TJ Pho- 
tons produced by radiation from an incoming or outgoing quark are called "prompt" ; an 
additional class of high-energy photons comprises those radiated from the incoming or 
outgoing lepton. In this letter, results are presented from a new inclusive measurement of 
isolated photon production in neutral current DIS. The data provide a test of perturbative 
QCD in a kinematic region with two hard scales: Q 2 , the exchanged photon virtuality, 
and Ej,, the transverse energy of the emitted photon. Compared to the previous ZEUS 
publication [7j, the kinematic reach extends to lower values of Q 2 and to higher values of 
Ej,. The statistical precision is also improved. 

Leading-logarithm parton-shower Monte Carlo (MC) and perturbative QCD predictions 
are compared to the measurements. The cross sections for isolated photon production in 
DIS have been calculated to order 0(a 3 ) by Gehrmann-De Ridder et al. (GGP) [SHU]. 
A calculation based on QED contributions to the parton distributions has been made by 
Martin et al. (MRST) [TJ]. 

2 Experimental set-up 

The measurements are based on a data sample corresponding to an integrated luminosity 
of 320 ± 8pb~\ taken between 2004 and 2007 with the ZEUS detector at HERA. The 
sample is a sum of 131 ± 3pb _1 of e + p data and 189 ± 5pb _1 of e~p data0 with centre- 
of-mass energy y/s = 318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [13J. Charged parti- 
cles were tracked in the central tracking detector (CTD) [H] and a silicon micro vertex 
detector (MVD) [TS] which operated in a magnetic field of 1.43 T provided by a thin su- 
perconducting solenoid. The high-resolution uranium-scintillator calorimeter (CAL) [16] 
consisted of three parts: the forward (FCAL), the barrel (BCAL) and the rear (RCAL) 

1 Hereafter 'electron' refers to both electrons and positrons unless otherwise specified. 
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calorimeters. The BCAL covers the pseudorapidity range -0.74 to 1.01 as seen from the 
nominal interaction point. The FCAL and RCAL extend the range to -3.5 to 4.0. The 
smallest subdivision of the CAL was called a cell. The barrel electromagnetic calorimeter 
(BEMC) cells had a pointing geometry aimed at the nominal interaction point, with a 
cross section approximately 5 x 20 cm 2 , with the finer granularity in the Z-direction 
This fine granularity allows the use of shower-shape distributions to distinguish isolated 
photons from the products of neutral meson decays such as tt° — > 77. 

A three-level trigger system was used to select events online [T7] by requiring well isolated 
electromagnetic deposits in the CAL. 

The luminosity was measured using the Bethe-Heitler reaction ep — > ejp by a luminosity 
detector which consisted of two independent systems: a lead-scintillator calorimeter [IS] 
and a magnetic spectrometer [T9] . 



3 Event selection and reconstruction 

Events were selected offline by requiring a scattered-electron candidate, identified using 
a neural network [20]. The candidates were required to have a polar angle in the range 
139.8° < 9 e < 171.9° in order to ensure that they were well measured in the RCAL. The 
impact point (X,Y) of the candidate on the surface of the RCAL was required to lie out- 
side the region (±15 cm, ±15 cm) centred on (0,0) to ensure well understood acceptance. 
The energy of the candidate, E' e , was required to be larger than 10 GeV. The kinematic 
quantities Q 2 and x were reconstructed from the scattered electron by means of the rela- 
tionships Q 2 = — (k — k') 2 and x = Q 2 / (2P ■ (k — k')) where k [k') is the four-momentum 
of the incoming (outgoing) lepton and P is the four-momentum of the incoming proton. 
The kinematic region 10 < Q 2 < 350 GeV 2 was selected. 

To reduce backgrounds from non-ep collisions, events were required to have a reconstructed 
vertex position, Z vtx , within the range |Z vtx | < 40 cm and to have 35 < 5 < 65 GeV, where 
5 = ^2 Ei(l — cos 8i); Ei is the energy of the ith CAL cell, 0j is its polar angle and the sum 

i 

runs over all cells [21] . At least one reconstructed track, well separated from the electron, 
was required, ensuring some hadronic activity which suppressed deeply virtual Compton 
scattering (DVCS) [22] to a negligible level. 

Photon candidates were identified as CAL energy-flow objects (EFOs) [23] for which at 
least 90% of the reconstructed energy was measured in the BEMC. EFOs with wider 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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electromagnetic showers than are typical of a single photon were accepted to allow eval- 
uation of backgrounds. The reconstructed transverse energy of the EFO, Ej,, was re- 
quired to lie within the range 4 < Ej, < 15 GeV and the pseudorapidity, 7/ 7 , had to satisfy 
—0.7 < rf 1 < 0.9. The upper limit on the reconstructed transverse energy was selected to 
ensure that the shower shapes from background and signal remained distinguishable. 

To reduce the background from photons and neutral mesons within jets, the EFO was 
required to be isolated from reconstructed tracks and hadronic activity. Isolation from 
tracks was initially achieved by demanding AR > 0.2, where AR = a/ (A(J)) 2 + (Ar/) 2 is 
the distance to the nearest reconstructed track with momentum greater than 250 MeV in 
the 7] — plane, where is the azimuthal angle. Jet reconstruction was performed on 
all EFOs in the event, including the electron and photon candidates, using the kx cluster 
algorithm [21] in the longitudinally invariant inclusive mode [25] with R parameter set to 
1.0. Further isolation was imposed by requiring that the photon-candidate EFO possessed 
at least 90% of the total energy of the jet of which it formed a part. 

Each event was required to contain both an electron and a photon candidate. The in- 
variant mass of the hadronic system, Wx, is then defined by W x = (P + k — k! — p 7 ) 2 , 
where p 7 is the four- vector of the outgoing photon. A total of 15699 events were selected; 
at this stage the sample was dominated by background events. The largest source of 
background was neutral current (NC) DIS events where a genuine electron candidate was 
found in the RCAL and neutral mesons, such as 7r° and rj, decaying to photons, produced 
a photon-candidate EFO in the BEMC. 

4 Theory 

Two theoretical predictions are compared to the measurements presented in this paper. In 
the approach of GGP [10], the contributions to the scattering cross section for ep — > e^yX 
are calculated at order a 3 in the electromagnetic coupling. One of these contributions 
comes from the radiation of a photon from the quark line (called QQ photons; Fig. Qja,b) 
and a second from the radiation from the lepton line (called LL photons; Fig. Hb,d). In 
addition to QQ and LL photons, the interference term between photon emission from 
the lepton and quark lines, called LQ photons by GGP, is evaluated. For the kinematic 
region considered here, where the outgoing photon is well separated from both outgoing 
electron and quark, the interference term gives only a 3% effect on the cross section. This 
effect is further reduced to ~ 1% when e + p and e~p data are combined as the LQ term 
changes sign when e~ is replaced by e + . The QQ contribution includes both wide-angle 
photon emission and the leading q — > q'j fragmentation term. GGP have chosen to use 
CTEQ6L leading-order parton distribution functions J2Sj. The factorisation scales used 
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are Q for QQ events and max(Q, fiF,min) f° r LL events where fiF,min = 1 GeV. Parton- 
to-hadron corrections were not made, in view of technical issues in relating 2—7-2 and 
2—7-3 topologies, following the advice of the GGP authors. We note that others have 
taken a different view [8] . A naive study indicated the likely effect to be a reduction after 
hadronisation in predicted inclusive cross-sections of order 15%. 

In the approach of MRST [T2"||2"T]. a partonic photon component of the proton, j p , is intro- 
duced as a consequence of including QED corrections in the parton distribution functions. 
This leads to ep interactions taking place via QED Compton scattering, 7 p e — > 7c A mea- 
surement of the isolated high-energy photon production cross section therefore provides 
a constraint on the photon density in the proton. The model includes the collinearly 
divergent LL contribution, which is enhanced relative to that of GGP by the DGLAP 
resummation due to the inclusion of QED Compton scattering. The QQ component is 
not included in the MRST model, in which the transverse momentum of the scattered 
electron is expected to balance approximately that of the isolated photon. In the analysis 
presented here, such a constraint was not imposed. The theoretical uncertainties in the 
models have been estimated by varying the factorisation scales by a factor two. 

Since the MRST cross sections include the LL contribution of GGP to a good approxi- 
mation, but exclude the QQ, an improved prediction can be constructed by summing the 
MRST cross section and the QQ cross section from GGP [2TPB"]. The theory uncertainties 
are of the same order as those of the individual QQ and LL components. 

5 Monte Carlo event simulation 

The MC program Pythia 6.416 [2H] was used to simulate prompt-photon emission for 
the study of the event-reconstruction efficiency In Pythia, this process is simulated as 
a DIS process with additional photon radiation from the quark line to account for QQ 
photons. Radiation from the lepton is not simulated in this Pythia sample. 

The LL photons radiated at large angles from the incoming or outgoing electron were 
simulated using the generator Djangoh 6 [30], an interface to the MC program Hera- 
cles 4.6.6 [21]; higher-order QCD effects were included using the colour dipole model of 
Ariadne 4.12 [32]. Hadronisation of the partonic final state was performed by Jetset 
7.4 [33]. The small LQ contribution was neglected. 

The NC DIS background was simulated using Djangoh 6, within the same framework 
as the LL events. This provided a realistic spectrum of mesons and overlapping clusters 
with well modelled kinematic distributions and hence was preferred to single-particle MC 
samples for backgrounds, such as were used in the previous ZEUS publication [7]. 
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The MC samples described above contained only events in which Wx was larger than 

5 GeV. Isolated photons can also be produced at values of Wx less than 5 GeV in 'elastic' 
and 'quasi-elastic' processes (ep — > epj) such as DVCS and Bethe-Heitler photon produc- 
tion. Such events were simulated using the GenDVCS [M] and Grape-Compton [35] 
generators. The contribution of these elastic processes was negligible after the selections 
described in Section [3j 

The generated MC events were passed through the ZEUS detector and trigger simulation 
programs based on Geant 3.21 [36J. They were reconstructed and analysed by the 
same programs as the data. In addition to the full-event simulations, MC samples of 
single particles (photons and neutral mesons) were generated and used to study the MC 
description of electromagnetic showering in the BEMC. 

6 Extraction of the photon signal 

The event sample selected according to the criteria in Section [3] was dominated by back- 
ground; thus the photon signal was extracted statistically following the approach used in 
previous ZEUS analyses [2HUI?]. 

The photon signal was extracted from the background using BEMC energy-cluster shapes. 
Two shape variables were considered: 

yi Ei\Zi~ Z e luster| 

• the variable (SZ) = - — - n ^ g . — , where Zj is the Z position of the centre of the zth 

i 

cell, Cluster is the centroid of the EFO cluster, w ce \\ is the width of the cell in the Z 
direction, Ei is the energy recorded in the cell and the sum runs over all BEMC cells 
in the EFO; 

• the ratio / max of the highest energy deposited in any one BEMC cell in the EFO to 
the total EFO BEMC energy. 

The distributions of (SZ) and / max (after the requirement (SZ) < 0.8) in the data and the 
MC are shown in Fig.[2J The MC LL and QQ distributions have been corrected in each 
two-dimensional (rj, Et) bin using factors derived from the difference between simulated 
and real DIS electron data. The (SZ) distribution exhibits a double-peaked structure 
with the first peak at ~ 0.1, associated with the signal, and a second peak at « 0.5, 
dominated by the tt° — > 77 background. The / max distribution shows a single peak at 
w 0.9 corresponding to the photon signal, and has a shoulder extending down to ~ 0.5, 
which is dominated by the hadronic background. 

The number of isolated-photon events contributing to Fig. [2] and in each cross-section bin 
was determined by a x 2 fit to the (SZ) distribution in the range < (SZ) < 0.8 using the 
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LL and QQ signal and background MC distributions as described in Section [51 By treating 
the LL and QQ photons separately, one automatically takes account of their differing 
hadronic activity, (resulting in significantly different acceptances) and their differing (77, 
Et) distributions, (resulting in different bin migrations due to finite measuring precision). 

In performing the fit, the LL contribution was kept constant at its MC-predicted value and 
the other components were varied. Of the 15699 events selected, 4164 ± 168 correspond 
to the extracted signal (LL and QQ). The scale factor resulting from the global fit for 
the QQ photons in Fig.[2]was 1.6; this factor was used for all the plots comparing MC 
to data. The fitted global scale factor for the hadronic background was 1.0. The signal 
fraction in the cross-section bins varied from 21% to 62%. In all cross-section bins, the 
X 2 /n.d.f. of the fits was 2.1 or smaller. 

For a given observable Y, the production cross section was determined using: 

da = A QQ ■ jV(7Q Q ) da^ 
dY C-AY dY ' 

where N(-jqq) is the number of QQ photons extracted from the fit, AY is the bin width, 
£ is the total integrated luminosity, cr^ is the predicted cross section for LL photons 
from Djangoh, and v4qq is the acceptance correction for QQ photons. The value of -4qq 
was calculated using Monte Carlo from the ratio of the number of events generated to 
those reconstructed in a given bin. It varied between 1.2 and 1.7 from bin to bin. 

The fits employed in this analysis were performed using (5Z) because of the larger dif- 
ference in shape between signal and background for this quantity. Fits in terms of the 
/max distributions were performed as a cross-check and gave similar results. As a further 
cross-check, an algorithm from the previous ZEUS publication [7], which selects wider 
electromagnetic clusters as photon candidates, was used. This proved to be more sen- 
sitive to the modelling of calorimeter backgrounds. In every case where a satisfactory 
fit was obtained, good agreement with the principal method was found. The corrections 
to the MC photon-signal energy-cluster shapes gave changes to the results within the 
statistical uncertainties and were not further considered [37] . 

7 Systematic uncertainties 

The following sources of systematic uncertainty were investigated |37j : 

• the energy scale of the electromagnetic calorimeter (EMC) was varied by its known 
scale uncertainty of ±2% causing variations in the measured cross sections of typically 
less than ±2%; 
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• the dependence on the modelling of the hadronic background by Ariadne was in- 
vestigated by varying the upper limit for the (SZ) fit in the range 0.6-1.0, giving 
variations that were typically ±5% but up to +12% and -14% in the most forward if 
and highest-x bins respectively. 

The following sources of systematic uncertainty were also investigated and found to be 
negligible compared to the statistical uncertainty [37] : 

• variation of the EMC energy-fraction cut for the photon candidate EFO by ±5%; 

• variation of the Z vtx cut by ±5 cm; 

• variation of the upper and lower cuts on 5 by ±3 GeV; 

• variation of the AR cut used for track isolation by ±0.1; 

• variation of the track-momentum cut used in calculating track isolation by ±100 MeV; 

• variation of the LL-signal component by ±5%. 

All the uncertainties listed above were added in quadrature to give separate positive and 
negative systematic uncertainties in each bin. The uncertainty of 2.6% on the luminos- 
ity measurement was not included in the differential cross sections but included in the 
integrated cross sections. 

8 Results 

The cross section for inclusive isolated photon production, ep —> cyX, was measured in 
the kinematic region defined by: 10 < Q 2 < 350 GeV 2 , W x > 5 GeV, E' e > 10 GeV, 
139.8° < 6 e < 171.8°, -0.7 < rf < 0.9 and 4 < < 15 GeV, with isolation such that 
at least 90% of the energy of the jet containing the photon belongs to the photon, where 
jets were formed according to the fcf algorithm with R parameter set 1.0. The measured 
integrated cross section is 

19.4 ± 0.7 (stat.)^; 2 , (syst.) pb, 
with an extracted contribution for QQ of 

12.2 ± 0.7 (stat.)tj; 2 (syst.) pb. 

The differential cross sections as functions of Ej,, ?7 7 , Q 2 and x are shown in Fig. [3] 
and given in Tables HHD It can be seen that the cross section decreases with increasing 
Erp, rf , Q 2 and x. The predictions for the sum of the expected LL contribution from 
Djangoh and a factor of approximately 1.6 times the expected QQ contribution from 
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Pythia agree well with the measurements, except for some differences at the lowest Q 2 
(and correspondingly lowest x). 

The theoretical predictions described in Section 2] are compared to the measurements in 
Fig. HJ The predictions from GGP describe the shape of the Ej, and 77 7 distributions 
well, but their central value typically lies 20% below the measured cross sections. The 
calculations fail to reproduce the shape in Q 2 \ a similar observation was made by HI [8]. 
As with the MC comparison, the measured cross section is larger than the theoretical 
prediction; this is also reflected in an excess of data over theory at low x. 

The MRST predictions mostly fall below the measured differential cross sections. How- 
ever, they lie close to the measurements at large values of Q 2 and x, for backward rp and 
for high values of Ej*, where the LL cross section is expected to be a substantial fraction 
of the total. Also included in Fig. |4] is the sum of MRST and QQ of GGP; it gives an 
improved description of the data over much of the range of the kinematic variables. 

Figure [5] shows the measured dajdrp compared to previous measurements from ZEUS [7] 
and HI [8] for the restricted range Q 2 > 35 GeV 2 and 5 < E^ < 10 GeV. The results are 
consistent but the uncertainty in the present measurement is smaller. 

9 Conclusions 

Inclusive isolated photon production has been measured in deep inelastic scattering using 
the ZEUS detector at HERA using an integrated luminosity of 320 pb _1 . Differential cross 
sections as functions of several kinematic variables are presented for 10 < Q 2 < 350 GeV 2 
and Wx > 5 GeV in the pseudorapidity range —0.7 < rp < 0.9 for photon transverse 
energies in the range 4 < Ej, < 15 GeV. The order a 3 predictions of Gehrmann-de 
Ridder et al. reproduce the shapes of the experimental results as functions of transverse 
energy and pseudorapidity, but are lower than the measurements at low Q 2 and low 
x. The predictions of Martin et al. mostly fall below the measured cross sections but 
are close in the kinematic regions where lepton emission is expected to be dominant. 
An improved description of the data is obtained by appropriately combining the two 
predictions, suggesting a need for further calculations to exploit the full potential of the 
measurements. 
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Ej, range 
(GeV) 


^(pbGeV- 1 ) 


4 6 
6 8 
8 10 
10 15 


4.87 ±0.28 (stat.) t° Q f 3 (syst.) 
2.40 ± 0.16 (stat.) IS?? (syst.) 
1.24 ±0.11 (stat.) (syst.) 
0.55 ± 0.04 (stat.) ±g;g| (syst.) 



Table 1: Measured differential cross-section -Mr, 



rf 1 range 




-0.7 -0.3 
-0.3 0.1 
0.1 0.5 
0.5 0.9 


17.4 ±0.9 (stat.) (syst.) 
13.0 ±0.8 (stat.) (syst.) 
10.7 ±0.9 (stat.) ±g;J (syst.) 
8.7 ±0.9 (stat.) ±1$ (syst.) 



Table 2: Measured differential cross-section -f^. 



Q 2 range 
(GeV 2 ) 


^(pbGeV- 2 ) 


10 20 
20 40 
40 80 
80 150 
150 350 


0.414 ± 0.035 (stat.) t tl (syst.) 

0.279 ± 0.020 (stat.) ±8:824 (syst.) 

0.115 ± 0.008 (stat.) IgjjJJ (syst.) 
0.050 ± 0.003 (stat.) Ztgloos (syst-) 
0.0088 ± 0.0009 (stat.) 1°;°°°^ (syst.) 



Table 3: Measured differential cross-section 4^. 



x range 




0.0002 0.001 
0.001 0.003 
0.003 0.01 
0.01 0.02 


5560 ± 380 (stat.) ±|° (syst.) 
3920 ± 230 (stat.) j^jjg (syst.) 

819 ± 58 (stat.) t% (syst.) 

103 ± 16 (stat.) t\l (syst.) 



Table 4: Measured differential cross-section 
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Figure 1: Lowest- order tree-level diagrams for isolated photon production in ep scat- 
tering. 
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Figure 2: Distributions of {5Z) and / max - The error bars represent the statistical 
uncertainties. The light shaded histogram shows a fit to the data of three components 
with fixed shapes as described in the text. The dark shaded histogram represents the QQ 
component of the fit, and the white histogram the LL component. The / max distribution 
is shown after requiring (5Z) < 0.8. 



21 



ZEUS 



> 

(5 
B 

Q. 
UJ 

e 

T3 



ZEUS 320 pb 1 
-LL + QQ(x1.6) MC 
■QQ(x1.6) MC 
■ LL MC 



(a) 



I I t LJ I I I 

10 12 14 

Ej(GeV) 



Q. 

*V 
e 

T3 




> 

O 
JD 

oj 
O 

T3 

D 
T3 



10 1 



(c) = 



10 



20 30 40 50 



100 



200 300 
Q 2 (GeV 2 ) 



J2 




S 6000 


X 




T3 

D 


5000 


T3 






4000 




3000 




2000 




1000 



io- 3 



(d) = 



Figure 3: Isolated photon differential cross sections in (a) E2,, (b) rp , (c) Q 2 and 
(d) x. The inner and outer error bars show, respectively, the statistical uncertainty and 
the statistical and systematic uncertainties added in quadrature. The solid histograms are 
the Monte Carlo predictions from the sum of QQ photons from Pythia normalised by 
a factor 1.6 plus Djangoh LL photons. The dashed (dotted) lines show the QQ (LL) 
contributions. 
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Figure 4: Data points as Fig. O Theoretical predictions from Gehrmann-De Ridder et 
al. and Martin et al. are shown with their associated uncertainties indicated by the shaded 
band and the hatched bands respectively. The dash-dotted line illustrates the combination 
MRST plus GGP: QQ. 
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Figure 5: Isolated photon differential cross-section compared to previous mea- 
surements at HERA with the additional kinematic restraints Q 2 > 35 GeV 2 and 
5 < Ej- < 10 GeV. The histograms show the different binnings used by ZEUS and 
HI. The symbols are mutually displaced for clarity. 
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